Berezinskii-Kosterlitz-Thouless transition in trapped quasi-two-dimensional Fermi gas 

near a Feshbach resonance 
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We study the superfluid transition in a quasi-two-dimensional Fermi gas with a magnetic field 
tuning through a Feshbach resonance. Using an effective two-dimensional Hamiltonian with renor- 
malized interaction between atoms and dressed molecules, we investigate the Berezinskii-Kosterlitz- 
Thouless transition temperature by studying the phase fluctuation effect. We also take into ac- 
count the trapping potential in the radial plane, and discuss the number and superfluid density 
distributions. These results can be compared to experimental outcomes for gases prepared in one- 
dimensional optical lattices. 
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I. INTRODUCTION 

The study of superfluidity and superconductivity in 
two-dimensional (2D) Fermi systems has attracted great 
attention in the past several decades, partly because of 
its close relationship to the problem of high- Tc super- 
conductors P, [1]. Recently, the experimental progress 
on creating quasi-low-dimensional atomic gases in opti- 
cal lattices [3I3,H,H01 and on atom chips provides 
us a possibility of realizing and studying superfluidity in a 
more controllable platform. In particular, with the aid of 
tuning an external magnetic field through a Feshbach res- 
onance, the interaction between fermionic atoms can be 
tuned continuously from the Bardeen-Cooper-Schrieffer 
(BCS) limit to the Bose-Einstein condensation (BEG) 
limit, so that the BGS-BEG crossover can be studied p^ . 

The BGS-BEG crossover in a 2D Fermi system has 
been discussed at both zero [ll[ and finite tempera- 
tures [13, where an effective 2D Hamiltonian with 
atom-atom interaction (called model 1) is employed. 
However, for strongly interacting fermionic atoms in a 
realistic quasi-2D geometry, due to the inevitable pop- 
ulation of many excited levels along the strongly con- 
fined transverse direction flSl, one needs to introduce 
a composite particle called dressed molecule to account 
for the population in these transverse excited levels, and 
write down a more complicated form for the effective 2D 
Hamiltonian (called model 2) with renormalized interac- 
tion between atoms (in the transverse ground level) and 
dressed molecules [15| . We have shown that the model 1 
and the model 2 Hamiltonians lead to qualitatively dis- 
tinct results in some cases .16!]. For instance, for quasi- 
2D fermions in a weak global harmonic trap, the model 

1 Hamiltonian predicts a constant value of the Thomas- 
Fermi cloud size at zero temperature from the BGS to the 
BEG limit. In contrast, the model 2 Hamiltonian predicts 
that the cloud size shrinks as one approaches the BEG 
limit [16], which is the correct trend as one should expect 
from the BGS-BEG crossover picture. 

In this paper, we extend the discussion of the model 

2 Hamiltonian to finite temperatures by taking into ac- 



count the fluctuation effects. In particular, we consider 
a trapped quasi-2D Fermi gas and study the superfluid 
transition temperature across a wide range of Feshbach 
resonance. Note that accordirig to the Mermin-Wagner- 
Hohenberg-Goleman theorem [13], phase fluctuations are 
dramatic at any finite temperatures such that there will 
be no ordinary long-range order in two dimensions. In 
this case, the superfluid transition is not accompanied 
by a Bose condensation, but instead characterized by 
a Berezinskii-Kosterlitz-Thouless (BKT) type, for which 
the vortex-antivortex pairs start to form at the transi- 
tion and a topological order is established [11] . Using the 
model 2 Hamiltonian, we show that in a uniform quasi- 
2D Fermi gas, the BKT transition temperature Tbkt in- 
creases from zero at the BGS limit and approaches mono- 
tonically to almost a constant value about 0.075£'f at the 
BEG limit, where Ep is the 2D Fermi energy. Further- 
more, we take into account the weak harmonic trap in the 
radial plane, and investigate the number and superfluid 
density distributions under the local density approxima- 
tion (LDA). As a characteristic feature of a BKT tran- 
sition, a finite jump of superfluid density is present in 
the middle of the trap, which takes a universal value of 
2rBKT/7r oM- 



II. FORMULATION 

A quasi-2D Fermi gas is typically prepared in exper- 
iments by means of a one-dimensional (ID) deep op- 
tical lattice along the axial (z) direction, where tun- 
neling between different lattice sites is completely sup- 
pressed. A ID lattice is described by the potential 
Vol ~ — Vocxp[— 2(x^ + y^)/w'^]cos{kzz), where is 
the wave vector of the laser beam, and w is the waist 
width satisfying w ^ This optical lattice poten- 

tial can be approximated around the minimal points by 
a strongly anisotropic pancake-shaped harmonic poten- 
tial with V « m^iUzz'^ +L0±r^)/2, where ujz = ^/Vok^/m, 

iLj± = 2y/VQ/{mw'^), and r — x"^ + y'^ is the radial dis- 
tance. The strong anisotropy of the trap {lOz 3> w±) 
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introduces two well separated energy scales, which al- 
lows us to first deal with the axial degrees of freedom 
by deriving an effective 2D Hamiltonian, while the radial 
degrees of freedom are left for later treatment with the 
LDA approximation. 

The effective 2D Hamiltonian describing a quasi-2D 
Fermi gas within a strongly axial confinement is derived 
in Ref. [T^- Under natural units {h — kg — 1), the 
Hamiltonian takes the form 
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where at^o- and dq are the annihilation operators for 
fcrmionic atoms and bosonic dressed molecules, respec- 
tively. The dressed molecule consists of a tight pair of 
atoms distributed in the transverse (axial) excited levels 
as well as the population in the closed Feshbach chan- 
nel [l^. In Eq. (1), eic = k^/(2m) is the dispersion 
relation for fermions with mass m and radial momen- 
tum k = {kx, ky), /i is the fermionic chemical potential, 
a =t, i labels the pseudo-spin, and is the quantization 
area. In the remainder of this manuscript, we choose lUz 
as the energy unit so that H, ii, and Ck = a^k^/2 become 
dimensionless, where = ^Jl/{mLOz) is the characteris- 
tic length scale for axial motion. 

The 2D effective bare parameters Ab (detuning), ajj 
(atom- molecule coupling rate), and Vb (background in- 
teraction) in Eq. [T] are connected with the 2D physical 



parameters Ap, 
relation 



and Vp through the renormalization 
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which is an identity for the variable x. The inverse of the 
physical effective interaction [T/g''*^(x)] are determined 
from the three-dimensional (3D) atomic scattering data 
through [l^ 



Ur, H — I - Sp{x) -f (7p{x) 

(3)" 

where Up = inabg/az, gp = i^coWUp/Uz, and jp = 
/ico(B — Bq)IlOz are 3D dimensionless parameters. Here, 
abg is the background scattering length, /ico is the differ- 
ence in magnetic moments between the open and closed 
collision channels, W is the resonance width, and Bq is 



the resonance point. The functions in Eq. ([3]) are [lE 
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where r(a::) is the gamma function. 

The partition function Z of the Hamiltonian Eq. ^ at 
temperature T = j3~^ can be expressed as an imaginary- 
time functional integral with action 



S 



dr 



E«k,a^rak,, 
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, (6) 



where a\^^„ — a\^^„[T) and rfq = dq(T) become imag- 
inary time dependent variables. By introducing the 
Hubbard-Stratonovich field &q(T), which couples to a^a^, 
and integrating out the fermions, we obtain Z = 
/ 6, d\d] exp(— ^cff), with the effective action given 
by 



|6qWP^ 
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Notice that the trace Tr is taken over momentum, imagi- 
nary time, and Nambu indices of the inverse Nambu ma- 
trix G^^, with 



G 



k,k' 
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(8) 

Here, Aq(r) = {ab/L)dq{T) + foq(T) is the field whose 
stationary value corresponds to the order parameter. 

We expand the action around a stationary and homo- 
geneous saddle point Ao = {ab/L)do + bo. The saddle 
point equation is derived from the saddle point action 
S'o[Ao] via the stationary conditions SSo/SAo = 0. This 
process leads to 
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(9) 

where i?k = \/ (ck — m)^ + Aq is the quasi-particle exci- 
tation spectrum. We have used the 2D renormalization 
relation Eq. Q in the derivation, which introduces the 
term l/(2ek + !)• 

In order to study the fluctuation around the saddle 
points of rfo, ^0, and Aq, we next write the field operators 



rfq(T) 

Aq(r) 



(Mo| + '5K(r)|)e'*<.M, 
i\bo\+S\bq{T)\)e^^-<^^\ 
(|AoK<5|Aq(r)|)e»M^) 



(10a) 
(10b) 
(10c) 
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which are all assumed to be slowly varying in both spa- 
tial and temporal coordinates. Here, S\dq\, <5|6q|, and 
S\Aq\ are magnitude fluctuations, while i^q, Tyq, and 6'q 
are phase fluctuations. The magnitude fluctuations cor- 
respond to gapped excitations, whose contribution can be 
neglected. However, the phase fluctuations correspond to 
gapless excitations, and it is their presence that breaks 
the ordinary long-range order in two dimensions. Notice 
that since the order parameter Aq is a linear combina- 
tion of dq and 5q, neglecting its magnitude fluctuation 
(5|Aq| is equivalent to assuming a same phase factor for 



dq and 6q. We thus have 



Mr) ^ Mr) EE 0q(r). 



(11) 



The phase fluctuation 6'q(r) around the saddle points 
is in general not a small value, so the phase factor in the 
order parameter cannot be directly expanded, but must 
be treated as a whole. This can be done by introducing 
a unitary transformation UO = exp[i6'q(r)/2]C' for the 
fields O = dq, 6q, and Aq. As a result of this trans- 
formation, and after a Fourier transformation back to 
the coordinate space, we can expand the effective action 
around its stationary point of Aq, leading to 



S'eff = ^o[Ao] + S'fluc[V0,5,6 



(12) 



where the fluctuation contribution S'duc contains terms 
only involving spatial or temporal gradients of the field 
6, which are assumed to be small. 

Expanding Sfiuc up to the quadratic order of W9 and 
(9r6' leads to 



S'fluc -^J'^drj d^r [iJdrO + K{drt 



-Pso(V0)2] , 
(13) 

where p^o represents the phase stiffness (or called the 
superfluid density) and takes the form 



PsO 



1 — tanh ^ ^ ^ 
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Here, Ub — A'^d[Vp^{2^)] ^/d^ denotes the dressed- 
molecular population. Other parameters in Eq. p3p are 
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By decomposing the phase fluctuation 6'(r, r) = 
Oyir, t) -|-0sw(r, t) into a vortex part 6^ and a vortex- free 
spin- wave part 0sw, one can separate the effective action 
S'cff, and hence the thermodynamic potential as 



(17) 
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FIG. 1: (Color online) Superfluid transition temperature 
Tbkt (solid lines) of a uniform quasi-2D Fermi gas for various 
values of Ef/oJz- Notice that Tbkt increases monotonically 
from the BCS (right) to the BEC (left) side of the resonance, 
and approaches to a Umiting value of ~ 0.075-Ef. In the BEC 
limit, Tbkt can also be estimated by considering a weakly in- 
teracting gas of composite bosons with 3D scattering length 
of 0.6aa, leading to a logarithmically varying result (dashed 
line) . The mean-field transition temperature Tmf (dashed line 
in the inset) is also shown as a comparison for Ef /ujz ~ 0.01. 



where fig, and flsw are the corresponding contribu- 
tions from the saddle point action 5*0, the vortex contri- 
bution from 9v, and the spin- wave contribution from 9sw, 
respectively. The spin-wave part can be integrated out to 
give ilsw = rX^fcli^ (l ~ 6"'=''/^), where c = y^plajj is 
the speed of the spin-wave. The number equation hence 
can be obtained via N = —dil/dn, leading to 



N 
1? 



2nb 
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-— tanh 

Sk V 2 



(18) 



with the spin- wave contribution rigw — —dflsw/dfi. By 
writing the number equation as above, we implicitly as- 
sume a vanishingly small vortex density, such that the 
vortex contribution to the number equation can be safely 
neglected. This condition is in general valid for temper- 
atures around or below the superfluid transition temper- 
ature f20'|. 

The is the superfluid density at the microscopic 
scale, and the vortex fluctuation tends to renormalize 
(decrease) the superfluid density when the physical scale 
I goes up. When the superfluid density is renomal- 
ized to zero with I —^ oo, the Berezinskii-Kosterlitz- 
Thouless (BKT) transition takes phase where the system 
changes from a superfluid phase to a normal phase. The 
renormalization-group (RG) flow depends on the vor- 
tex core energy Ec, which can be estimated by Ec = 
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FIG. 2: (Color online) Two-body binding energy \E\,\ (solid 
lines) as a function of a^/as for various values of Ef/uIz- No- 
tice that \Eb\ can exceed Ef (dotted) even on the high-field 
side of the resonance point, hence pushes the system into the 
BEC limit as indicated by the transition temperature shown 
in Fig. [U 



0.787r/5so 1211 ■ The RG equations for the superfluid den- 
sity ps and the fugacity y — e^^^" are given by [2^ 



dK-\l) 
di 

dyjl) 
dl 



= 4^V(0 + O(y3), (19) 
= [2-i:K{l)]y{l) + 0{y^), (20) 



where K{1) = ps{l)/T. The initial conditions are ps{l — 
0) = Pso a-iid y{l = 0) = exjp(—TTpsQ/2T) . The fixed point 
{y{oo), K{oo)) — (0, 2/7r) of these RG equations corre- 
sponds to the critical condition for the BKT transition, 
where the vortex- antivortex pairs start to dissociate and 
break superfluidity. Thus, the equation for the critical 
temperature Tc — Tbkt is determined as [l^l 



BKT 



(21) 



where pf — liuii^oo Ps{l) is the renormalized superfluid 
density. This equation must be solved together with the 
saddle-point equation ([5]), the number equation (jlSp . and 
the RG equations (fTO|) and ^U\i to determine p, Aq, and 
Tbkt as functions of the magnetic field detuning. As a 
comparison, the mean-field transition temperature Tmf 
can be calculated by solving only Eqs. ([9]) and (flSj) with 
Ao = 0. Next, we will consider some specific systems and 
discuss the corresponding numerical results. 



III. NUMERICAL RESULTS 

We first ignore the trapping potential in the radial 
plane, and consider a uniform quasi-2D Fermi gas. A 



typical set of results of Tbkt for various values of ujz 
are shown in Fig. [TJ Here, we consider the case of ^Li 
and use the parameters ai,g — — 1405ao, W = 300 G, 
and pco — '^Pb, where ag and ps are Bohr radius and 
Bohr magneton, respectively. We also notice that when 
plotted as functions of the inverse 3D scattering length 
Qz/as, the results for '^'^K are very close to those for ^Li, 
indicating the universal behavior around the resonance 
point. 

Notice that the superfluid transition temperature 
Tbkt increases continuously from the BCS to the BEC 
side of the Feshbach resonance, and saturates to a lim- 
iting value of Tbkt ~ 0.075Ef. This limiting value of 
Tbkt is expected for psQ sa n/{Am), which corresponds 
to the BEC limit of a quasi-2D Fermi gas where paired 
fermions behave like weakly interacting bosons with num- 
ber density n/2 and mass 2m. In the BEC limit, we can 
also consider the particles as composite bosons with ef- 
fective three-dimensional scattering length of 0.6as [23| . 
For weakly interacting bosons, more accurate results are 
available for the BKT transition temperature based on a 
combination of the quantum Monte Carlo simulation [23| 
and the renormalization group approach 25], which pre- 
dicts that Tbkt goes down logarithmically with the scat- 
tering length. This more accurate result in the BEC limit 
is also shown in Fig.[Tl which is slightly deviant from our 
fermionic calculation. 

We note also from Fig. [T] that the BEC limit value of 
Tbkt can even be reached on the BCS side of the Fesh- 
bach resonance (with < 0). In fact, since a two-body 
bound state is always present in quasi-low dimensions at 
arbitrary detunings, the BEC limit of a quasi-2D Fermi 
gas can be realized as long as the condition \Et\ :S> Ep 
is satisfied, regardless of the sign of the 3D scattering 
length Og. This property is in clear contrast to the 3D 
case, where a bound state is only present with positive 
as, hence the BEC limit can only be realized on the low- 
field side of a Feshbach resonance. Here, the two-body 
binding energy \Ei,\ is determined by solving the Hamil- 
tonian ([1]) for two particles, which gives 



[VfiE,)]' 



1 



4TTa'. 



■H-Eb). 



(22) 



The corresponding results for \Eb\ are plotted in Fig. 
Notice that for the case with Ep/uj^ — 0.005, the binding 
energy \Eb\ « 2.QEp for Oz/ag ^ —2.0, leading to a BKT 
transition temperature of the BEC limit value Tbkt ~ 
0.075£;f, as shown in Fig.[TJ 

Another feature of Fig.[T]is that the transition temper- 
ature Tbkt increases with the axial trapping frequency 
LOz at a fixed detuning. This trend can also be easily un- 
derstood by analyzing the two-body binding energy 
In fact, at a given detuning, |i?f,| increases with ojz as 
shown in Fig. [21 hence pushes the quasi-2D system fur- 
ther to the BEC hmit. 

The results of Tbkt should also be compared with the 
outcome of an effective 2D Hamiltonian with renormal- 
ized atom-atom interaction (model 1), which is discussed 
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FIG. 3: (Color online) Superfluid transition temperature 
Tbkt (solid lines) plotted as functions of the two-body bind- 
ing energy \Ei,\. The results are compared with the out- 
come of an effective 2D Hamiltonian with renormalized atom- 
atom interaction (model 1, dotted line), which is discussed in 
Refs. 



in Refs. [13. fl3j. In Fig. [31 we show Tbkt calculated us- 
ing both models as functions of the binding energy \Eb\. 
One prediction of the model 1 is that the many-body 
physics (such as the BKT transition temperature) takes 
a universal behavior, which only depends on the two- 
body binding energy \Eb\ in the unit of Ep [HI, d, El] ■ 
As shown in Fig. [31 this is clearly not the case for the 
results from the model 2, where the transition tempera- 
ture Tbkt also depends on the other energy scale such as 
the transverse trapping frequency w^. Notice that both 
of the models predict roughly the same limiting value of 
Tbkt for large values of \Et,\ (in the BEC limit), since 
in that limit the system behaves like weakly interacting 
bosons, and Tbkt gets very insensitive to the interaction 
between these composite bosons (logarithmic dependence 
as mentioned above). 

After discussing the case of a uniform quasi-2D Fermi 
gas, we next consider the trapping potential in the ra- 
dial plane V±{r) — mu!±r'^/2. Under the local density 
approximation (LDA) , we consider a position dependent 
chemical potential /i(r) — fio — V±i'i^)- Here, /io is the 
chemical potential at the trap center, which must be de- 
termined by fixing the total particle number in the trap 
N = 271 Jp°° n{r)rdr. Using this technique, we calculate 
the in-trap number and superfluid density distributions 
in the radial plane, which are shown in Fig. [H From 
the left to the right, results are sequentially presented 
for the BCS side of the resonance (a^/cs = —1), uni- 
tarity {uz/as — 0), and the BEC side of the resonance 
{az/cLa = 1). For each case, temperature is varied within 
two orders of magnitude, showing that superfluidity is 
absent at high temperatures and starts to build from the 
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FIG. 4: (Color online) In-trap number (solid) and superfluid 
(shaded area) density distributions at (a) T — O.Oluz, (b) T — 
O.OOlojz, and (c) T = O.OOOlojz. For each temperature, three 
typical detunings are considered with Oz/oa = — 1 (BCS, left), 
dz/cis = (unitarity, middle), and a^/aa = 1 (BEC, right), 
respectively. Notice the finite jump of superfluid density in 
the middle of the trap at intermediate temperatures. The 
plots are made for ^Li, while the results for ''"K are similar. 
Parameters used in this plot are lo^ = 2-n x. 62 kHz, lu±_ = 
2-K X 10 Hz, N = 10000, and Rbcs = '^/2lJ7JujI{Nf/'^az is 
the zero-temperature Thomas-Fermi cloud size of a quasi- 2D 
ideal Fermi gas with particle number A'". 



trap center as decreasing T. At low enough temperatures, 
superfluidity nearly extends to the whole trap, and the 
number density profile approaches the zero temperature 
results [l6l|. Furthermore, notice that there is a finite 
jump of superfluid density present in the middle of the 
trap at intermediate temperatures, as shown in panels of 
Fig. ^p. This discontinuity in signatures the position 
where the superfluid transition takes place, and takes a 
universal value of 5p^ = 2T /it characterizing a phase 
transition in the BKT universality class J^]. 



IV. CONCLUSION 

In summary, we have studied the BKT type of su- 
perfluid transition of a quasi-2D Fermi gas based on an 
effective 2D Hamiltonian with renormalized interaction 
between atoms and dressed molecules. The finite tem- 
perature effect has been taken into account by incorpo- 
rating phase fluctuations over the saddle point solution. 
Using the effective Hamiltonian that we derived before 
with the explicit parameters [l^, we establish the BKT 
transition temperature as a function of the 3D atomic 
scattering length, making it possible to directly compare 
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the result with the experimental measurements. We also 
compare the predictions from this effective Hamiltonian 
(model 2) with the results from model 1 (where the 2D 
Hamiltonian is described by atom-atom interaction with 
an effective scattering length ) , and they differ sig- 

nificantly in certain parameter regions. In particular, the 
universal behavior predicted by the mode 1 (many-body 
physics depends only on the two-body binding energy) 
is not the case for the model 2, where the BKT tran- 
sition temperature also has explicit dependence on the 
transverse trapping frequency lOz. This difference can be 
quantitatively tested by future experiments. Under the 



local density approximation, we have also calculated the 
in-trap number and superfluid density distributions at 
finite temperatures, which can be compared with exper- 
imental results in a harmonic trap. 
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